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Abstract In a simplified model we study how close can the result of a fast non-equilibrium process, like hadronization of 
quark matter in an expanding fireball, come to hadronic composition in equilibrium. In a chemical approach of in-medium 
quark fusion (A + B —* C) we find that this depends on the time-integrated rate per volume multiplied by the conserved 
number of constituents. Ideal hadron gas equilibrium constrained by conservation laws, the fugacity parametrization, as well 
as linear and non-linear coalescence approaches are recognized as different approximations to this chemistry. It is shown that 
color confinement requires a dependence of the hadronization cross section on quark density. 
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Since strangeness enhancement had been suggested as a possible signature of quark-gluon plasma formation in 
heavy-ion collisions, there is a particular attention to the question of hadronic and quark level equilibrium in these 
reactions 0-^]. Recent claims of describing data at AGS as well as at CERN SPS by constrained hadronic equilibrium 
(i.e. using a common volume, temperature and only chemical potentials of conserved quantities for the computation 
of hadron numbers) "quite" satisfactorily [^|-|9| led to some astounishment in the community: the total reaction time 
' at SPS energy is namely too short for achieving hadronic equilibrium by hadronic conversion processes alone. Pushing 
aside the faux pais, that this "satisfactory" equilibrium description agrees with experimental data for a few ratios 
only within a factor of two while the general agreement is on the 30% level, one faces the question, how can this 
C^l | simple idea work even so well. The time argument above would exclude this. 
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Meanwhile there is a number of non-equilibrium approaches to this problem ranging from the use of phcnomeno- 
logical parameters describing the deviation on the hadronic level (fugacities) [10,[llj through linear jl2| and nonlinear 
(constrained) coalescence models |l3|,|l4]], and mixed quark-hadron chemistry 16 1 up to very detailed microscopic sim- 



ulations of the hadronization (parton cascades fll7| i20|l VENUS |21| , HIJING |22 



B3j, RQMD E4E5L UrQMD [E6| M 



string model [p0|-p3| , chromodielectric model |34| ) . The basic physical concepts of these approaches oppose hadronic 



equilibrium; it might be reached in course of time, but for some components it is actually not. For recent transport 
simulations at SPS energies see p5[ . There are also recent doubts whether chemical equilibrium can be reached at 
I 1 RHIC ||§. 

states with quark level processes instead of hadronic ones in a fast expanding fireball. In order to concentrate on this 
quantifying effort of the time evolution of a complicated system we shall deal with a simple fusion process on the 
constituent level: 



A + B < — ► C. 

An example is given by the fusion of an s quark and a ud diquark to A particle. In this case baryon charge and 
strangeness conservation constrains the equilibrium. 



In such a process the energy-impulse balance is assumed to be satisfied by the surrounding medium, we are interested 
in the time evolution of the number of final state products Nc{t). The equilibrium number of precursors NA, eq and 
Ns,eq - in a quark matter scenario - will eventually be taken as zero. The initial state is assumed to be far from 
this equilibrium, it is oversaturated by A and B. The simplest is to consider no C particles at all in the beginning. 
As the reaction proceeds there will be C type particles at the end (or equilbrium will never be reached), they are 
hadrons. The cross section for such a reaction - also as a simplification - can be taken as a constant parameter, but 
temperature and reaction volume must be considered as time dependent quantities. We shall point out, however, 
that the assumption of constant cross sections actually contradicts to the confinement principle in quark matter 
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scenarios, so in that case the hadronization cross section has to be quark density dependent. This extra dependency, 
not derivable from perturbative or lattice QCD so far, will be referred as "dynamical confinement" in this paper. 



In the followings - after a brief presentation of general chemical equations - we shall consider the simple model 
of hadronization due to coalescence (fusion in the presence of a medium). Comparison of different limiting cases 
(allowing or excluding the precursors A and B in equilibrium) and estimates of the characteristic deviation from 
hadronic equilibrium with and without dynamical confinement will close this study. 



In the chemical approach to this process all particle numbers are time dependent. The entropy grows during a 
spontaneous approach towards equilibrium. A general reaction can be characterized by stochiometric coefficients a[ r ^ 
in a fashion arranged to zero: 

i 

(r) 

for each reaction r. In this arrangement the a\ -s can be negative as well as positive. An important quantity with 
respect to approaching chemical equilibrium is the activity of a reaction, 

i 

It describes the (un)balance of gain and loss terms. Rate equations summarize contributions from each reaction 
channel, 

N^Y.a^R^il-e^). (2) 

(r) 

Here the factors depend on numbers of other particles, cross sections and volume. Due to known (e.g. adiabatic) 
expansion a relation between volume V and temperature T is given. 



In simple scenarios of expansion VT 3 / 2 = const, (massive ideal gas expands adiabatically, no back effect from 
chemistry is taken into account). Now 

Ni = N i>eq e^ T (3) 

and 

Nw = Vdi(^y 2 e- m '' T . (4) 

In simple scenarios like this the time dependence of particle numbers, Ni(t), is entirely due to that of the 
(Hi(t) — m,i)/T{t). The equilibrium ratio, 

R = A^.eqA^eq = d A d B / m A m B \ 3/2 3/2 ( mc -mA-m B )/T /r-v 

° q JVc,eq ' d c \ m c J [ ' 

is nevertheless constant for adiabatic expansion (I^T 3 / 2 =const.) and composite final products (mc = m A + ms). 



All rate equations drive the system towards chemical equilibrium. The total entropy production can be obtained 
from 

TdS + VidNi = dE+ P dV = Q (6) 
for adiabatic flow of an ideal ("dry") fluid. One obtains 

S = - VJ ^N, = VJ BW S W (> m - l) > o. (7) 
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In chemical equilibrium — for each reaction, i.e. 

E^ = °- ( g ) 

i 

There are as many linear constraints for the chemical potentials (ii as many reactions. If this is less than the number 
of particle sorts, the equilibrium state is undetermined. If this is greater than the number of particle sorts, the 
reactions cannot be all independent: there are conserved quantities (e.g. from Ni = N2 it follows that B = N\ — N2 is 
conserved). In this latter case all ^ can be expressed by a few numbers corresponding to these conserved quantities. 
This is the constrained equilibrium. 



Let us denote the conserved quantities by Q a (a = baryon, strangeness, etc.). The corresponding charge of the i-th 
particle sort let be qi >a . We have 

E<7i,JVi = (9) 

Substituting the rate equations (||) we get 

i r 

Since conservation is fulfilled even out of equilibrium, we conclude that 

'£4 r) q i , a = Q^=0, (11) 

i 

i.e. there is no charge unbalance in any of the reactions. This fact allows us to redefine the chemical potentials as 

Mi =y^gi,aMa + M»- ( 12 ) 

a 

The first contribution, depending on only as many "chemical potentials" as conserved charges we have, does not 
change the activities: 

^EE^^+E^I- (13) 

i a i 

Since the first part is zero due to eq. (pT[) , chemical equilibrium (vanishing activities) means 

E«i r) f = 0- (14) 

i 

For a well-determined or overdetermined system the only equilibrium solution is given by fa = 0. In this (constrained) 
equilibrium state the original chemical potentials are a linear combination of a few /i a -s, associated to conserved 
charges, 



:,e 9 = E 9i > a ^ a - ( 15 ) 



Whether this chemical equilibrium state can be reached in a finite time depends on the rates RS r ' and activities 
x^ r \ In this respect whether a quantity is conserved or not is determined by the set of rate equations. This way it 
happens that not only real conserved charges but also the number of each constituent quark may be conserved in a 
sudden hadronization process, like ALCOR assumes this. 
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For ideal massive particles (Boltzmann distribution) 

Mi = Tlog— ^. (16) 

The activities hence contain ratios of equilibrium numbers in each reaction (mass action law, detailed balance) . For a 
simplified expansion picture these ratios are (nearly) constant. The factors contain precursor numbers, averaged 
cross sections (fluxes) and the reaction volume. In an expanding system its time-dependence is far from trivial. 

For a general set of reactions and particle sorts the time evolution can only be followed numerically. In order to 
understand, however, possible deviations from equilibrium in a final state - even after infinitely long time - we consider 
an oversimplified system in the followings. 



In the case of only one reaction channel, A + B < — ► C, we denote the equilibrium ratio by R eq (cf. eq.(|5j)). Since 
ola = olb = 1 and ap = —1 for this very reaction, we have Na = Nb = —Nc- From this two conservation laws follow, 

N a + N c =Na(0) = N , (17) 

N B + Nc = N B (0) =N . (18) 
For considering an evolution towards equilibrium the specific rate, 

(crv) 

^=V' (19) 

is a useful quantity. It occurs in the rate equation 

N c = v (N A N B - R eq N c ) • (20) 
Substituting the conservation laws (|l^JT^) we arrive at 

Nc = v {{N - Nc){N - N c ) - R eq N c ) = v{N c - N+){N C - JV_). (21) 
This rate is vanishing for the particle numbers Nc — N± with 



N± = R eq /2 + N ± ^JRIJA + R eq N . (22) 

Here iV_ is a stable, N + is an unstable equilibrium point. For a non-relativistic, adiabatic expansion 

VT 3/2 = const. (23) 

and hence R eq = const, for constituent hadrons (those with mass mc = rriA + tub)- A realistic limit is R eq = when 
no A or B particles are allowed in equilibrium. In particular this is the case for quark confinement. 



The simple rate equation (21) can be solved analytically giving 

rNc(oo) ^ . th 

J Nc{0) (N C - N + )(N C - N-) ~J t0 



i/dt = r. (24) 



Here to is the initial time, when the fireball contains no hadrons, Nc{0) — 0. t^ is the break-up time, beyond which 
no chemical reactions (even no particle collisions) occur. In some cases it is large compared to characteristic times 
and can be approximated by infinity. Now using the identity 

(25) 



(N c - N+)(N C - N-) N+-N-\N C -N + N c - N- 
we arrive at 
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log(iV c - N + ) ~ \og(N c - AL)| ^ = {N+ _ 7V_)r. (26) 

This leads to the solution 

Nc ^ = N + N - N+-KN. (27) 

with 

K = aq>{-(N + -N_)r). (28) 

Here N— is the constrained equilibrium value, achieved for K = 0. There are two alternative forms of this relation 
derived by using N + N- = Nq. The first one is the coalescence form, 

N c (oo) = iV 2 _i_A_ = C(K,N Q )N* = C(K,N )Na(0)N b (0), (29) 

the second one is the equilibrium form, 

1 -K 
1 - K(N-/N ) 



N c (oo) = iV_ T ^ rAT /Ar ^ 2 = lc N_ = lc N c , eq . (30) 



In particular Nc(oo) = N- in the K = (infinite rate) case. On the other hand in the case of confinement R eq = 0, 
from which N+ = iV_ = A^o follows. 

This agrees with the ALCOR assumption: after long enough time all quarks A and B are converted into hadrons 
C. This means A^_ — Nq, satisfied for R eq — 0, i.e. when the quark equilibrium numbers are zero (confinement in 
equilibrium). Then ALCOR ends up with the equilibrium number. (Note: ALCOR does not take into account the 
gain/loss balance as if it were infinite activity in favor of hadrons - against quarks. This feature is equivalent to 
putting R eq = 0.) 

Without quark confinement, but (N + — N-) small (0(R eq )), i.e. when precursors are suppressed in the equilibrium 
state, we consider a negligible, but not zero R cq . In the limit of R oq — > the precursors A and B are suppressed in 
equilibrium, according to the static confinement principle. This is, however, insufficient for hadronizing all quarks in 
an expanding matter. We show this in the followings. 

We expand K in terms of (A^ + — A^_) which is in the same order of magnitude. We get 

K = cxp(-(A r + - N-)r) = 1 -r{N+ - N_) (31) 
and the final hadron number becomes 



Eq. ( p2| ) - valid if the precursors are suppressed in chemical equilibrium (and this is all we know from lattice QCD) 
- has different readings. For r small it comes close to the linear coalescence equation 

7V c (oo) = rN% = rN A (0)N B (0). (33) 

(Note that Na(0) = Nb(0) = No.) ALCOR uses extra modification factors 6a and bs by assuming the modified 
coalescence rate 

Nc(oo) = rb A b B Nl (34) 
To start with it seems to be rather similar to the coalescence picture. Due to conservation, however, ALCOR obtains 

AT c (oo) = rb A b B Nl = —^N* (35) 
1 + rN 
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and determines the factors hi accordingly (depending on No). It leads exactly to the constrained equilibrium result 
in the limit rNo — > oo: 

N c (oo)=N . (36) 

So, reversing the argument, we conclude that for quark matter hadronization r has to diverge in order to reflect 
confinement in the end stage. 



Finally the chemistry approach gives the interpretation of the fugacity due to 

A c (oc) = Wo j^Jvo = lcN °- (37) 

In fact the difference between these approaches is more pronounced by inspecting chemical potentials instead of final 
numbers. In the discussed simple case the chemical potential of the (final state) hadron is given by 

tta Na (oo) 

^=log-^. (38) 
Substituting our result for Nc (oo) it splits to contributions with different physical meanings: 

~~T~ = l0g N^7 q + l0g ^v" = T + l0g7c ' (39) 

The conserved quantities (charges) carried by the hadron C are those involved in the identities respected by the rate 
equations, N A (t) + N c (t) = N A (0) = N and N B (t) + N c (t) = N B (0) = N . In this case q c ,i = 1 and q c ,2 = 1, 
respectively. The fugacity parameter (as used by Rafelski et al.) is the remaining part of the formula, in equilibrium 
it is 7c = 1. 



Finally we note that the ALCOR approach (constrained coalescence) leads to equilibrium result only if a single 
reaction channel is considered, like in the case discussed in the present paper. For several, competing channels the 
relative hadron ratios are different for ALCOR and constrained equilibrium. This best can be seen by inspecting the 
corresponding fa chemical potentials to the pre resonance decay hadron numbers as a function of strangeness in the 
ALCOR simulation. Slight deviations from the perfect straight lines, which would mark the constrained equilibrium 
in the ideal case, can be observed p7| . 



In case of confinement for particles A and B the equilibrium ratio R eq vanishes and A_ has to be replaced by Nq 
in the above formula. Exactly this way are usually the chemical potentials associated to conserved charges, \i\ and 
fj,2, determined. The fugacity factor 



1-K 

1C ~ 1 -K(N-/N ) 2 ' ( ' 



in only static confining scenarios reduces to 



1 + rN 



This must be viewed as an indication for the necessity of considering dynamic confinement in hadronization. 



There is no difference between ALCOR and constrained equilibrium in this simple model if confinement is im- 
plemented: both lead to jc = 1 eventually. The linear coalescence picture on the other hand assumes rN -C 1, 
corresponding to a rather impressive deviation from equilibrium. Therefore it is interesting to estimate the factor rNo 
for different heavy-ion experiments. 



We continue by estimating the integrated specific rate r. Its value depends not only on cross section(s) - first 
assumed to be constant and later density dependent - but also on the expansion and cooling of the reaction zone. 



G 



We consider here two limiting cases: that of the one dimensional expansion with a rather relativistic equation of 
state (ideal gas with particles having negligible masses) and that of three dimensional expansion of a massive, non- 
relativistic ideal Boltzmann gas. The first case is thought to be rather characteristic for RHIC and LHC, the latter 
rather for SPS and AGS energies. 



During a one-dimensional expansion the volume grows linearly in time, V = itR^t, and the (thermal) average of 
relative velocities is near to the light speed (v) = 1. In this case we obtain 

a dt (ctvq) t br 

The time-integrated specific rate grows with the break-up time tb r monotonically, it can eventually be infinite, if there 
is no thermal break-up of the expanding system. In that case 7c = 1, equilibrium is achieved. For a three-dimensional 
expansion one gets V — Vb(i/io) 3 - Furthermore a massive, non-relativistic ideal gas satisfies VT 3 ^ 2 — VqTq^ 2 during 
an adiabatic expansion. The average relative velocity is temperature dependent like (v) ~ y/T/m ~ 1/t, with some 
reduced mass of the fusing pair m. This together leads to the estimate 

^ r '"4 d 4- (43) 



Evaluating the integral we arrive at 




It is interesting to note that r is finite even without break-up, i.e. in the tb r — > 00 limit. This also means that the non- 
relativistic, 3-dimensional spherical expansion scenario with constant hadronization cross section cannot comply with 
quark confinement. In this case a has to be time dependent at least as a ~ i 3 ~ V ~ 1/n, i.e. the hadronization cross 
section has to be inversely proportional to the density of quarks. Exactly this assumption was made in transchemistry 

M. 



In orther to simplify the further discussion we relate the thermal average of the rate at the beginning of the 
expansion to estimates of the collision time due to mean free path, 

to = 7-A-. (45) 
{avo )pa 



with density po = 2Nq/Vq. This leads to 

(<7V ) 1 



V 2N t c 

The quantity describing the deviation from constrained equilibrium for hadrons then is given as follows: 



(46) 



riV = ^log^ (47) 

Zt c to 



rNo=£-U-(±) ! (48) 



for one-dimensional relativistic flow and 

^ 3 

for three-dimcsional non-relativistic flow of massive constituents. The general dependence of the fugacity on the 
estimated collision time is in both cases 

1 + p{t c /t ) 
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just the constant j3 differs. For the realtivistic, one-dimensional flow case it is 



= — (50) 



and for three-dimensional spherical flow of massive particles it is given by 



P= 7-T3 

From these formuli one realizes that for tbr — to, i.e. in case of immediate break-up the chemical equilibrium cannot 
be approached, no hadron C will be formed and 7c remains zero. The opposite limit, tbr — > 00 is less uniform: while 
in the relativistic one-dimensional scenario equilibrium is always reached if enough time is given (7(7(00) = 1), this is 
not so for a stronger, more spherical expansion. In the latter case 7c tends to a finite value, 

^ (oo) = TT^Ao' (52) 

which is in general less than one. It means that chemical equilibrium cannot be really reached even in infinite time. 
(An analysis of proton - deuteron mixture with a similar conclusion was given in Ref. [p8|.) For early break-up, 
tbr — to *C to, both scenarios lead to 

0=-^. (53) 

tbr — to 



In quark matter hadronization scenarios eq.(p2|) contradicts to the quark confinement principle. In order to reach 
7c (00) = 1 a dynamical confinement mechanism, causing the quark fusion cross section in medium to be (at least) 
inversely proportional to the quark density, has to be taken into account besides the equilibrium suppression of quark 
number (referred in this paper as "static confinement" ) . 



Let us assume that A and B are quark level (colored) objects, and the fusion cross section scales as 

a = ao (L 2 ^o V +e (54) 
°\V N A + N B J V ' 

with some extra power e. The solution of the rate equation in the R cq = case (both static and dynamic confinement) 
becomes 

N c = N Q (l - (1 - S) 1 ^ , (55) 

with 

< f/t - ■ 



For e = — 1 (without dynamical confinement) one gets back the result discussed so far. In the limiting case e = 
(minimal dynamical confinement) one obtains 

to 



N c = iV 1 - ^ . (57) 



Fig.l plots the dependence of the non-equilibrium fugacity factor 7c on the scaled break-up time tbr /to for these 
scenarii: one and three dimensional scaling expansion, in the three-dimensional case both for static and dynamic con- 
finement (color density dependent quark fusion cross section). The uppermost curve, showing the fastest equilibration 
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belongs to a 3-dimensional expansion scenario with a cross section proportional to the inverse of the quark density 
(cf. eq.([57|)). The lowest curve, saturating well below the equilibrium value, corresponds to the same expansion with 
static confinement only. The middle curve shows the static RHIC scenario, reaching equilibrium (vanishing quark 
number) eventually, but quite slowly. For all calculations we used the ratio to /t c = 3 for the sake of simplicity. (The 
smaller to, the higher density, therefore the shorter inter-collision time t c initially.) 

Equilibration 

1 I 1 1 1 1 1 1 1 1 1 



0.9 - SPS dynamic 




1 2 3 4 5 6 7 8 9 10 

tor 't 

FIG. 1. Degree of the chemical equilibrium reached for different break-up times tbr/io and initial inter-collison times t c /to 
for SPS and RHIC expansion scenarios with static and dynamic confinement. 



Finally we compare characteristic values for SPS and RHIC. In the SPS scenarios we assume t = 1 fm/c as 
the starting time of the expansion. The collision time changes somewhat due to different temperatures (and effective 
quark masses) presumed at the beginning of hadronization. We consider a value characteristic for quark-gluon plasma, 
t c » 0.33 fm/c. We take over the ratio to/t c = 3 for the RHIC scenario, too. The break-up time is the most unknown 
factor, therefore the equilibration has been plotted as function of ibr/^o- Fortunately, once the stauration value of 
7c has been approached, there is no strong dependence on it. The 1 — e _1 ~ 0.63 fraction of equilibrium value is 
reached after 2 fm/c for the dynamic SPS, after 3.5 fm/c for the static RHIC scenario. In the static SPS scenario it 
will be never reached, jc saturates near 0.3 already after 2 fm/c. This last scenario, quark coalescence with density 
independent probability, can be excluded on the basis of quark confinement. 



In conclusion 1) the results of the statistical model can be interpreted not only by assuming hadronic equilibrium, 
but as well by quark - hadron mixture equilibrium with quarks eventually suppressed due to confinement. This 
interpretation avoids the collison time problem: the hadrons need not transmutate into each other directly. 2) The 
ALCOR and constrained hadronic equilibrium leads to the same result in this simple scenario due to too many 
constraints on conserved quantities. Quark coalescence, which takes conservation into account, therefore mimics 
hadrons in chemical equilibrium. This feature certainly does not not generalize for several competing reaction channels, 
but the overall picture can be similar. ALCOR actually assumes more conserved quantities, than usual: quark and 







antiquark numbers in any flavor are conserved separately due to the assumption of rapid hadronization. 3) jc ~ 0.3 is 
typical for CERN SPS without considering dynamical confinement, even for an arbitrary long hadronization process. 
This leads to a factor two to three characteristic error on equilibrium model results. The experimental evidence, 
reaching at least 70-80% of equilibrium ratios for almost all hadrons in a short time in the order of several fm/c, 
excludes this possibility. It also underlines the quark matter hadronization scenario with dynamical confinement 
effect, as transchemistry describes it. 4) An expansion scenario typical for RHIC and LHC energies in principle 
allows for reaching the equilibrium even with static confinement, but the process is rather slow and an early break-up 
may prevent from hadronizing all quarks. Dynamical confinement has to be in work for one-dimensional expansion 
scenarios, too. 



One is tempted to speculate about the origin of dynamical confinement, looking for a cause of the density dependence 
of the hadronization cross section. Since this dependence makes the cross section stronger with diminishing quark 
density, it cannot have perturbative origin; on the contrary it has to do with formation of strings or color ropes. As 
the color density drops, strings, follow-ups the chromoelectric flux lines, become longer before ending on a partner 
color. The volume, in which possible partner to form a hadron occur, hence increases, so does the hadronization cross 
section, a ~ £ 3 . The characteristic string length on the other hand depends on the color density as ~ n C o/or • Such a 
mechanism would lead to a = cro7i C oior(0)/n C oior. 
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